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I.  Introduction 


The  origins  of  this  position  paper  lie  in  discussions  about  the  proper  form  of  the  compliance 
matrix  and  its  inverse,  the  modulus  matrix,  appropriate  for  Terfenol-D.  With  very  little 
additional  work,  we  can  discuss  the  full  piezomagnetic  matrix  which  includes  not  only  the 
compliance  matrix  but  also  the  piezomagnetic  constants  and  the  permeability  constants.  The 
general  matrix  has  45  independent  constants,  however  these  can  be  reduced  to  5  independent 
constants  by  applying  symmetry  appropriate  to  today’s  high-power  Terfenol-D  drivers.  We  also 
briefly  discuss  methods  of  measuring  these  constants. 

II.  Table  of  symbols 

d  effective  piezomagnetic  constant  (A/m) 

e  effective  piezomagnetic  constant  (m/A) 

g  effective  piezomagnetic  constant  (T) 

h  effective  piezomagnetic  constant  (T_1) 

sH  elastic  compliance,  H  =  constant  (m2/N) 

sB  elastic  compliance,  B  =  constant  (m2/N) 

cH  elastic  stiffness,  H  =  constant  (N/m2) 

cB  elastic  stiffness,  B  =  constant  (N/m2) 

H  magnetic  field  strength  (A/m) 

B  magnetic  flux  density  (T) 

fd1  permeability,  stress  =  constant  (H/m) 

jus  permeability,  strain  =  constant  (H/m) 

vT  reluctivity,  stress  =  constant  (m/H) 

Vs  reluctivity,  strain  =  constant  (m/H) 

8,  S  strain 

T,  T  stress  (N/m2) 

E  Young’s  modulus  (N/m2) 

III.  General  Elastomagnetic  Matrix 


In  tensor  form,  the  four  basic  linear  relationships  for  the  piezomagnetic  materials  become: 
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In  these  relationships:  (1)  the  shear  strain  components  are  defined  by  %  =  ( l/2)(du/dxj  +  du/dxi), 
i  4  j,  where  the  m;’s  are  the  components  of  the  displacement  vector,  and  (2)  the  shear  stress 
components  are  not  compressed  may  not  equal  ty,  allowing  for  body  rotations). [1] 

It  is  an  accepted  procedure  to  compress  the  2nd  rank  strain  and  stress  tensors  to  6  dimension 
vectors  and  the  4th  rank  tensors  into  36  component  2nd  rank  tensors.  Employing  the  Voigt  strain 
components  defined  by  Sij  =  du/dxj  +  du/dxi  =  2 and  the  stress  components  given  by  Ty  =  Ty  + 
tji,  the  first  general  elastomagnetic  matrix  listed  above  can  be  written: [2,3] 
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The  remaining  three  elastomagnetic  relationships  follow  in  a  similar  manner,  but  they  will  not  be 
considered  here. 


Note  that  in  terms  of  the  true  tensor  notation,  this  matrix  becomes:  [1] 
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This  matrix  has  three  components: 

1.  The  6x6  s  submatrix  is  the  compliance  matrix;  it  has  21  independent  components. 

2.  The  3x6  d  submatrix  is  the  piezomagnetic  matrix;  it  has  18  independent  components. 

3.  The  3x3  ju  submatrix  is  the  permeability  matrix;  it  has  6  independent  components. 

IV.  Effect  of  Symmetry 

If  we  consider  an  isotropic  material  polarized  with  a  magnetic  field  and  uniaxial  stress 
applied  along  the  z  axis,  the  matrix  simplifies  greatly.  Recall  that  SS,  ST,  SB  and  SH  are  small 
quantities  in  the  elastomagnetic  formulation.  If  larger  values  of  stress  and  field  are  used,  as  may 
be  the  case  in  high  power  sound  projectors,  this  relationship  yields  average  values  that  may  be  a 
function  of  stress  and  field  amplitude.  For  this  case,  the  first  elastomagnetic  matrix  reduces  to: 
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Terfenol,  on  the  other  hand,  is  not  isotropic,  structurally  or  magnetically.  Structurally  the 
[112]  direction  is  the  preferred  axis  of  growth  with  twinned  platelets  of  Tbi_xDyxFe2  separated 
by  eutectic  sheets.  In  the  large  Terfenol  samples,  the  platelets  are  not  parallel  throughout  the 
sample  and,  to  the  lowest  approximation,  can  be  considered  randomly  oriented  in  the  sample 
plane.  Magnetically  the  magnetostriction  constants  are  not  isotropic,  in  fact,  Am  »  Am,  which 
makes  the  piezomagnetic  d  constants  very  anisotropic.  The  approximation  which  most  closely 
resembles  the  large  diameter  Terfenol-D  rods  in  the  transducer  configuration  consists  of  isotropy 
in  the  plane  perpendicular  to  an  applied  field  and  stress  H  and  T  along  the  growth  axis.  (To  date, 
it  is  impossible  to  prepare  high  quality  large  diameter  [112],  [111],  or  [110]  single  crystals.) 
Fortunately,  for  this  case  of  magnetic  field  and  stress  applied  along  the  [112]  axis  and  assuming 
isotropy  in  the  plane,  the  correct  matrix  is  identical  to  the  fully  isotropic  case  shown  above. [1] 

In  the  tensor  notation,  this  can  be  rewritten: 
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where  we  have  set  du  =  -  ^25  =  0,  see  Ref.  1. 

V.  Uniaxial  Transduction  Excitation 


If  in  addition  to  the  bias  stress  and  field  along  the  z  axis,  the  a.c.  fields  and  stresses  are 
applied  along  the  same  (z)  axis,  only  the  following  equations  remain: 


SSxx  =  s13H8Tzz  +  d318Hz 
dSyy  —  Sl^ST^  +  d3]SHz 
8S zz  -  s33h8Tzz  +  d338Hz 
8Syz  =  0 

8Sv  =  0 

<5^  =  0 
SBX  =  0 
8By=  0 

8BZ  =  d338Ta  +  fx33  §HZ 

At  this  stage  in  the  development  of  Terfenol,  it  is  not  possible  to  incorporate  more  terms  and 
characterize  the  material  better.  Large  diameter  Terfenol-D,  as  prepared  today  and  utilized  in 
currently  proposed  high  power  sonar  systems,  can  be  fully  characterized  by  the  measurement  of 
the  following  five  coefficients: 

1.  s33H=1/E, 

2.  sj3h  -  -vIE,  where  v  is  Poisson’s  ratio, 

3.  d33t 

4.  d3i, 

5.  JU33T, 

and  all  five  “constants”  are  functions  of  both  H  and  T. 
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VI.  Recommended  Measurement  Techniques 


Assuming  stress  and  field  are  applied  only  in  the  z  direction  and  that  the  plane  perpendicular 
to  the  stress  is  isotropic  in  the  plane,  Terfenol-D  can  be  characterized  in  a  straightforward 
manner:  Strain  gages  are  cemented  both  along  the  direction  of  applied  stress  and  field  and  in  the 
perpendicular  direction.  One  large  strain  gage  or  a  number  of  smaller  gages  can  be  used  in  order 
to  obtain  a  good  average  of  the  constants.  (Strain  gages  of  course  measure  only  the  surface  and 
it  must  be  assumed  that  the  material  does  not  change  radially.)  One  possibility  is  to  apply  a 
single  strain  gage  around  the  circumference  of  the  sample  or  average  a  number  of  smaller  gages. 
An  optical  method  or  linear  variable  transformer  method  (LVDT)  can  also  be  used  to  obtain  an 
average  over  the  entire  sample,  at  least  in  the  z  direction. 

The  applied  field  and  stress  can  be  measured  by  a  number  of  conventional  methods.  One 
method  is  to  employ  a  large  face  electromagnet  whose  pole  faces  are  modified  to  apply 
compressive  stresses  of  sufficient  magnitude.  A  load  cell  measures  the  applied  stress.  B  can  be 
measured  using  a  pickup  coil  wound  around  the  sample;  the  field  can  be  measured  by  a  Hall 
probe.  Alternatively,  dual  concentric  coils  can  be  would  around  the  sample.  By  measuring  the 
flux  through  both  of  the  coils,  both  B  and  H  can  be  calculated. 

1.  H  is  held  constant  and  T is  varied.  The  compliances  S3311  and  snH  and  Poisson’s  ratio  v 
can  be  obtained  from  the  parallel  and  perpendicular  strain  gages  respectively,  c/33  can  be 
obtained  from  this  measurement  if  B  is  measured  during  the  change  in  stress. 

2.  T  is  held  constant  and  H  is  varied.  The  piezomagnetic  constants  and  d\ 3,  are 
obtained  from  the  parallel  and  perpendicular  strain  gages  respectively.  The 
permeability  //33r  is  obtained  by  measuring  B. 

These  measurement  techniques  are  also  valid  for  the  lower  hysteresis  TbxDyyHozFe2 
quaternary  samples  that  are  prepared  using  the  same  processes  as  the  Terfenol-D  samples. 

VII.  Numerical  Examples 

As  discussed  earlier,  the  “constants”  in  the 
elastomagnetic  matrices  are  not  really  constants.  Fig.  1 
shows  the  variation  of  the  d  constant  and  ju  relative  (=  jul 
jUq)  as  a  function  of  H  for  a  large  diameter  Terfenol-D  rod. 

The  d  “constant”  varies  by  almost  a  factor  of  10  as  H  varies 
and  //  relative  varies  by  a  factor  of  2.5.  Fig.  2  shows  the 
compliance  s33H  as  a  function  of  T  at  a  field  of  147  kA/m 
(-1850  Oe).  s33H  varies  by  a  factor  of  1.8.  For  these 
numerical  examples  we  must  pick  a  field  and  a  stress  and 
find  the  appropriate  parameters.  For  our  field  H,  we  choose 
the  maximum  d33  point,  151  kA/m  (-1900  Oe),  d33  =  9.1  x 

10'9  m/A.  The  relative  ju33  at  the  same  point  is  3.8.  For  a 
compressive  stress  of  56  MPa,  s33H  is  3.2  x  10'11  m2/N.  We 
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use  the  values  of  Poisson’s  ratio  v=  0.45  and  the 
value  of  J33  /  J31  =  -2.0  given  by  Claeyssen[4],  This 
gives: 

1.  533"=  3.2  x  lO-Hm2^ 

2.  snH  =  -1.44x  lO-Hn^/N, 

3.  ^33  =  9.1  x  10'9  m/A, 

4.  J31  =  -4.6  x  10‘9  m/A,  and 

5.  p 337  =  3.8  piQ  =  4.8  x  10'6  T-m/A. 

Assuming  that  SHZ  =  36  kA/m  (-450  Oe)  and  ST  -  0, 
we  have: 

SSxx  =  =  £?3i  5Hz  =  -6.91  x  10‘5  m, 

SSZZ  =  ^33  <3?/z  =  3.26  x  10'4  m,  and 
SBZ  =  SHZ  —  0. 172  T 


Fig.  2.  Compliance  vs.  stress  for  a  large, 
Bridgman,  Terfenol-D  sample  at  147  kA/m. 


Obviously,  if  there  is  no  magnetoelastic  coupling,  which  corresponds  to  all  of  the  d’ s  being  zero, 
all  of  the  SS’s  are  also  zero. 

VIII.  Summary 

Large  diameter  Terfenol-D  samples  with  the  stress  and  field  applied  in  the  z  direction  and 
isotropic  properties  in  the  x-y  plane  can  be  described  by  only  five  parameters.  These  parameters 
can  be  measured  with  a  pickup  coil  and  strain  gages  mounted  parallel  and  perpendicular  to  the  z 
axis. 


If  the  preparation  methods  improve  in  the  future  such  that  high  power  transducer  material 
can  be  made  similar  to  the  FSZM  materials  or  [111]  and  [110]  single  crystals  can  be  produced, 
then  the  symmetry  in  the  x-y  plane  is  broken  and  a  larger  number  of  constants  must  be  used  to 
characterize  the  materials. 
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